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Resolving the nuclear dust distribution of the Seyfert 2 
galaxy NGC 3081 
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ABSTRACT 

We report far-infrared (FIR) imaging of the Seyfert 2 galaxy NGC 3081 in the range 70- 
500 [xm, obtained with an unprecedented angular resolution, using the Herschel Space 
Observatory instruments PACS and SPIRE. The 11 kpc (~70") diameter star- forming 
ring of the galaxy appears resolved up to 250 |i.m. We extracted infrared (1.6-500 |J.m) 
nuclear fluxes, that is active nucleus-dominated fluxes, and fitted them with clumpy 
torus models, which successfully reproduce the FIR emission with small torus sizes. 
Adding the FIR data to the near- and mid-infrared spectral energy distribution (SED) 
results in a torus radial extent of Ro=4±f pc, as well as in a flat radial distribution of 
the clouds (i.e. the q parameter). At wavelengths beyond 200 |a.m, cold dust emission 
at T=28±l K from the circumnuclear star-forming ring of 2.3 kpc ('^15") in diameter 
starts making a contribution to the nuclear emission. The dust in the outer parts of 
the galaxy is heated by the interstellar radiation field (19±3 K). 

Key words: galaxies: active - galaxies: nuclei - galaxies: imaging - galaxies: indi- 
vidual (NGC 3081). 



1 INTRODUCTION 

The infrared (IR) SED of active galactic nuclei (AGN) 
serves as a sensitive probe of both the dust and the 
sources that are heating it. Dust grains absorb opti- 
cal and ultraviolet photons from the AGN and from 
stars and re-radiate them in the infrared (IR) range. 
The IRAS and ISO satellites revealed that Seyfert galax- 
ies are strong FIR and mid-infrared (MIR) em itters 
ilRodriguez Espinosa et al.l 119871 : ISpinoglio et al.l [l995|) and 
that this emission is thermal, and a combin ation of a warm, 
a cold, and a very cold dust components I Ra.dovich et al.l 



1 1993 : IPerez Garcia fc Rodriguez Espinosall200ll ) . The warm 
component is produced by dust heated by either the AGN or 



E-inail:C. Ramos@sheffield.ac.uk 



circumnuclear starbursts (dust at 120-170 K), the cold dust 
is heated by stars in the galaxy disk (30-70 K), and the very 
cold dust is heated by the general interstellar radiation field 
(15-25 K). With the advent of the Herschel Space Obser- 
vatorjij (|Pilbratt et al.ll2010l ) it is now possible to map the 
FIR emission of nearby Seyferts at higher angular resolu- 
tions than those previously achieved with the Spitzer Space 
Telescope between 70 and 160 (xm. This, together with the 
unprecedented sensitivities that Herschel offers up to 500 
\x.m allow to probe their dust distributions at different tem- 
peratures. In this letter we present new Herschel imaging 



^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 
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data of the galaxy NGC 3081, which is part of a guaranteed 
program of FIR imaging observations of Seyfert galaxies. 
The main goal is to characterise their IR SEDs, by deter- 
mining the fractional contributions of the warm, cold, and 
very cold dust components. 

The ga laxy NGC 3081 ha rbours a Seyfert 2 nucleus, al- 
though iMoran et al.l ((20001) reported a spectacular Type-1 
optical spectrum i n polarised light. Th e galaxy is at a dis- 
tance of 32.5 Mpc (jButa &i Purcelllll998l ). which corresponds 
to a spatial scale of 158 pc arcsec"^. This early-type barred 
spiral is forming stars in a series of nested ringlike features: 
a nuclear 2.3 kpc diameter ring (hereafter rl), an inner ring 
of 11 kpc (r2), an o uter ring of 26.9 kpc, and a pseudoring of 
33.1 kpc di ameter JButa fc Purcelllll998l : ISuta et al.ll2004l : 
iBvrd et al.ll2006l '). 

2 OBSERVATIONS 

FIR maps of NGC 3081 were obtained with the PACS and 
SPIRE instruments of the Herschel Space Observatory. The 
data are part of the guaranteed time proposal "Herschel 
imaging photometry of nearby Seyfert galaxies: testing the 
coexistence of AGN and starburst activity and the nature 
of the dusty torus" (PI: M. Sanchez-Portal). 

The PACS observations were carried out using the 
"mini-map" mode, consisting of two concatenated 3' scan 
line maps, at 70° and 110° (in array coordinates). This re- 
sults in a map with a highly homogeneous exposure within 
the central 1' area. The PACS beams at 70, 100, and 
160 ^m are 5.6", 6.8", and 11.3" full-width half maxi- 
mum (FWHM) respectively. With the SPIRE photometer, 
the three available bands were observed simultaneously us- 
ing the "small map" mode, whose area for scientific use is 
around 5'x5'. The FWHM beam sizes at 250, 350, and 500 
\im are 18.1", 25.2", and 36.9" respectively. 

We carried out the data reduction with the Herschel In- 
teractive Processing Environment (HIPE) v6. 0.1951. For the 
PACS instrument, we deemed the extended source version 
of the standard PhotProject reduction script as adequate, 
given the small angular size of the gala xy. We used the EM 
v5 photometer response calibration files (|Miiller et al.ll2011f ). 
For SPIRE, we applied the standard small map script with 
the 'naive' scan mapper task, using th e calibration database 
v6.1. Colour correct ions (for PACS, seelPoglitsch et al.ll2010l : 
please refer to the ISPIRE OmI |2011| for the SPIRE ones) 
are small for blackbodies at the expected t emperatures (e.g. 
IPerez Garcia fc Rodriguez Espinosal [200l|) and have been 
neglected. More details on the observations and data pro- 
cessing are given in Sanchez- Portal et al. (in preparation). 

The FIR maps of NGC 3081 are shown in Figure [T] 
The star-forming ring r2 is clearly resolved in the three 
PACS images, as well as in the 250 ^im SPIRE map. There 
is a brighter region in the western side of the ring de- 
tected at 70, 100, 160, and 250 \im that does not have 
optical/near-infrared ( NIR) counterpar t in the HST imag es 
iButa fc Purcelllll998l : rButa et al.|[20o3 : lBvrd et al.ll2006l '). 

To study the IR nuclear emission of NGC 3081, we ob- 
tained and compiled unresolved fluxes, i.e. either the emis- 
sion of a point spread function (PSF) component fitted to 
the data, or the emission contained in an aperture diameter 
equals to the FWHM of the PSF in each band. In the case 
of the FIR, we used GALFIT 2D fitting (|Peng et al.l 120021 ') 
to obtain the unresolved fluxes. For the PACS images we 
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0.047 
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0.35 


0.055 
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T-ReCS 


18.3 


0.56 


0.083 


231 


25 


PACS 


70 


5.6 


0.83 
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30 


PACS 


100 


6.8 


1.01 


575 


30 


PACS 


160 


11.3 


1.7 


60 


30 


SPIRE 


250 


18.1 


2.7 


488 


50 


SPIRE 


350 


25.2 


3.7 


86 


50 


SPIRE 


500 


36.9 


5.4 


19 


60 



Table 1. IR nuclear fluxes employed in the fit of NGC 3081 
with clumpy torus models. Errors have been obtained by adding 
quadratically the photometric and PSF subtraction uncertainties. 



fitted a PSF component, which we identified with the nu- 
clear flux, a Sersic profile of Re ~ 7" — 18", and a fainter 
and larger Sersic component of Re ^ 45" — 60", where Re 
is the half-light radius given by GALFIT. For the SPIRE 
data we only fitted a PSF component and a Sersic profile of 
Re ~ 75" - 90". We also tried this simpler PSF + Sersic fit 
with the PACS images, but the three-component model re- 
sults in smaller residuals and values of the reduced x^ (0.002 
for PACS and 0.02 for SPIRE). All the fitted Sersic compo- 
nents have indices between 0.4 and 0.8 (i.e. disk-like) and 
do not reproduce the r2 ring, which is a residual of the fits. 
The PSF in put functions are the empirical on es from PACS 
and SPIRE (|Lut j|2010l : Isibthorpe et al.ll201ll ). 

In addition, we compiled the highest angular resolu- 
tion NIR and MIR data from the literature to construct 
the nuclear SED. Subarcsecond resolution MIR images of 
NGC 3081 (0.30" at 8.74 ^m and 0.56" at 18.3 ^m) were 
obtained using the camera/spectrograph T-ReCS on the 
Gemi ni-South Telescope. The u nresolved T-ReCS fluxes 
from iRamos Almeida et al.l ((20091) are reported in Table [TJ 
together with an additional nuclear flux at 13.04 ^.m from 
VISI R on the VLT with similar resolution as the T-ReCS 
data ((Gandhi et al.l 12009). In the NIR, we use the nuclear 
flux obtained from the NICMOS 1.6 (im image reported in 
iQuillen et all ((200 ll ). 

In Table [T] we report the unresolved NIR, MIR, and 
FIR fluxes, their uncertainties, and the angular resolution 
at each wavelength. The FIR errors are the result of adding 
quadratically the photometric accur acies (5% for PAC S at 
70 and 100 ^.m and 10% at 160 >j .m; iMiiller et"ai]|2010l and 
7% for SPIRE: ISPIRE OM||2011| ') and the PSF flux deter- 
mination uncertainties. The latter are the dominant source 
of error and account for the variations in the PSF fluxes 
associated to the GALFIT fitting in each band. 

3 NUCLEAR SED MODELLING 

In a recent series of papers ((Ramos Almeida et al.l (2009|; 



2011: lAlonso-Herrero et al.ll201ll ') we fitted the nuclear NIR 
and MIR emis sion of Sey fert galaxies with the clumpy torus 
models of Ne nkova et al.l ( 2008) and were able to nicely con- 
strain the torus parameters, including the torus radial ex- 
tent, Y = Ro/Rd, where Ro and Rd are the outer and inner 
radius of the toroidal distribution of clumps respectively. 
The inner radius is defined by the dust sublimation temper- 
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Figure 1. Herschel PACS 70, 100, and 160 i^m images (top) and SPIRE 250, 350, and 500 |J.m maps (bottom). North is up and East is 
to the left. The r2 ring is resolved up to 250 |xni. 



ature (T^j ~ 1500 K). However, it is not clear if the lack of 
FIR high angular resolution data, which probes cooler dust, 
might bias the fits to smaller torus sizes. In the context of 
the Nenkova models, clouds are heated by the AGN radi- 
ation (directly illuminated) and by other clouds (indirectly 
illuminated) and each clump contains a range of tempera- 
tures itself. The temperature of the dust within the torus 
scales with the square root of the distance to the sublima- 
tion radius. Thus, a distribution of clumps with Y ~ 100 will 
include dust at temperatures ranging from Td to ~150 K for 
the directly illuminated clouds, and down to a few Kelvin 
for the shadowed and more distant ones. To test how the ad- 
dition of FIR data affects the fits described above, here we 
use Herschel PACS and SPIRE nuclear fluxes combined with 
NIR and MIR data of the galaxy NGC 3081 and fitted them 
with clumpy torus models. In general, clumpy torus models 
appear to reproduce better the IR emission of nearby AGN 
than smooth torus m odels (e.g. lAlonso-Herrero et al.ll2003l : 
iMullanev et al.ll2011f ) although there is not a general con- 
sensus on the dust distribution yet. Indeed, in forthcoming 
publications based on guaranteed time Herschel observations 
of AGN we plan t o use the smooth torus mod e ls des cribed in 
iFritz et all (|2006l ) and lHatziminaoglou et all (|2008l ) to com- 
pare with the results obtained with clumpy torus models. 

We constructed the nuclear 1.6-500 ^.m SED of NGC 
3081 (see Figure [2| using the nuclear fiuxes reported in Ta- 
ble [T] The unresolved NIR and MIR components correspond 
to a physical region of <85 pc in diameter and thus, are likely 
be dominated by emission from the AGN dusty torus. The 
nuclear FIR fluxes of NGC 3081, on the other hand, come 
from regions with sizes of ~1 kpc in the case of PACS and 



between 2.7 and 5.4 kpc for SPIRE. These regions are much 
larger than the physical scales responsible for the NIR and 
MIR unresolved emission. Indeed, the SPIRE nuclear fluxes 
include emission from the inner ring of 2.3 kpc diameter 
(rl), and consequently, we consider them as upper limits 
in the fit (see Figure [21l. The PACS fluxes exclude rl, but 
may include other sources of nuclear emission apart from the 
torus. This contamination might affect the SED shape and 
consequently, the resulting torus parameters. However, here 
we work under the assumption that the torus is the domi- 
nant source of u nresolved emission up to 160 (im. Although 
iDeo et al.l l|2009l) showed that, in general, the AGN contin- 
uum of Seyfert 2 galaxies drops rapidly beyond 20 (xm, NGC 
3081 is one of the galaxies in their sample with the smallest 
starburst to AGN ratios at 30 \ira. Indeed, the Spitzer IRS 
flux measurement that they reported for this galaxy, 1.09 
Jy, nicely matches our fitted torus models (see below). 

The clumpy dusty torus models of iNenkova et al.l 

(120081 ) are characterised by six parameters, which are de- 
scribed in Table [2] Here we use an interpolated version 
of the Nenkova models to fit the IR nuclear fluxes re- 
ported in Table [l] (considering the SPIRE fluxes as up- 
per limits) using our Bayesian inference tool BayesClumpy 
(JAsensio Ramos &: Ramos Almeidal |2009| ) and the uniform 
priors described in Tabled The result of the SED fitting are 
the posterior distributions of the model parameters (Figure 
[Sjl, but we can translate these results into a single SED. In 
Figure [2] we plot the model that better fits the IR data, i.e. 
the maximum-a-posteriori (MAP) model, and the model de- 
scribed by the medians of the six posteriors resulting from 
the fit (see Table [2]). The nuclear NIR, MIR, and FIR emis- 
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Figure 2. Rest-frame IR SED of NGC 3081 (dots). Solid and dasiied lines are the "best fit" to the data (MAP) and the model described 
by the median of the posteriors respectively. The shaded region indicates the range of models compatible with the observations at the 
2(T level. The SPIRE fluxes are set as upper limits in the fit. 
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Figure 3. Normalised marginal posteriors (NMP) resulting from 
the fit of NGC 3081. Solid and dashed vertical lines represent the 
modes and medians of each NMP respectively, and dotted vertical 
lines indicate the 68% confidence level for each parameter around 
the median. 



sion of NGC 3081, at least up to 160 |i.in, is successfully re- 
produced by a relatively broad clumpy torus (cr=57°) with 
an average number of clouds No=5 along the radial equato- 
rial direction, and with a close to edge-on incl ination i=71°. 
In Ramos Almeida et al. (2011; hereafter I I RAl 11) we fit- 



ted the NIR/MIR nuclear SED of NGC 3081 using the same 
models as here and obtained a, No, and i values which agree 
with those reported here. On the other hand, the torus radial 
extent, Y, and the index of the clouds radial distribution, q, 
change significantly when adding the FIR data. The sensi- 
tivity of the SED to q for small values of Y is highly reduced. 
This is because for small tori the SED shape does not change 
noticeably either when the clumps are distributed along the 
whole extent of the torus {q=0) or highly concentrated in its 
inner part (q=2-3). Here we obtain g=0.2, which is charac- 
teristic of a flat clo ud dist ribution, and a torus radial extent 
y=84, whereas inl lRAlll the resulting values for NG C 3081 
were 5=2.3 and y=22. However, while in | IRAIII we im- 
possed V=[5, 30] as a prior, in this work we use Y=[5, 100] 
to take into account cold dust within larger scales. We did 
the test of fitting only the NIR and MIR data using the same 
priors shown in Table [2] and obtained y=35, with the other 
parame ters resulting in similar values as those reported in 
I Ira 111 . We also performed the fit by considering both the 
PACS and SPIRE fluxes as upper limits, and the results are 
practically the same as those from the NIR/MIR flt. Sum- 
marising, including the FIR data in the flt of NGC 3081 
results in a relatively large torus radial extent and flattens 
the clouds distribution. 

Using the median value of Ay (in the following we will 
refer to median values in order to give uncertainties at the 
68% confidence level, which are defined around the median) 
reported in Table [2] we can derive the column density using 



.9 X lO-"' 
4.1±?:ixl0^ 



the Galactic dust -to-gas ratio {Nj^ 

iBohlin et aLlll978l ). This gives ]<ijf^ = 4.1±i:6 XIU"" cm 
which is compatible with the value derived from ASCA 
X-ray obse rvations of NGC 308 1 (N^"™"" = 6.3 ± 0.4 x 
10^^ cm~^: iLevenson et al.ir2009l ). The AGN bolometric lu- 
minosity can be obtained from the vertical shift applied to 



Resolving the nuclear dust distribution of NGC 3081 5 



Parameter 



Abbreviation Interval 



Fitting results 
Median Mode 



Width of the angular distribution of clouds 

Radial extent of the torus {Ro/Rd) 

Number of clouds along the radial equatorial direction 

Power-lawr index of the radial density profile 

Inclination angle of the torus 

Optical depth per single cloud 



o- 


[15°, 75 


Y 


[5, 100] 


No 


[1, 15] 


q 


[0,3] 


i 


[0°, 90° 


TV 


[5, 150] 



48°±{2„ 

76±1| 
8±| 

o-4±g:t 

54°±"° 

25 

59±if 



57° 

84 

5 

0.2 

71° 

72 



Ay produced by the torus along the line-of-sight (LOS) 



j^LOS 



214:±ll'' mag 



162 mag 



Table 2. Clumpy model parameters and Ay^ derived from them. Columns 1 and 2 give the parameter description and abbreviation 
used in the text. Column 3 indicates the input ranges considered for the fit (i.e., the uniform priors). Finally, columns 4 and 5 list the 
medians and modes of the posterior distributions shown in Figure \3\ 



the models to fit the data: L^^j 



jxlO^^ erg s~ 



= 2.1±J;| 

and can be directly compared with the 2-10 keV intrinsic 
luminosity derived from the ASCA data, after applying a 
bolometric correction factor of 20 (Elvis etal,! 1994) ; L^; — 
1.0±0.2xl0''^ erg s'^ The differenc e betw een L^jf^ and 
L^j is smaller than the one found bvl lRAllI for NGC 3081. 
As explained before, the outer size of the toroidal distribu- 
tion of clouds is defined as Ro = YRd, where Kd scales with 
Lft". Thus, Ro = 0.4 Y (L«/«/10«)0-5 pc = 4±? pc, 
which includes dust at very different temperatures (from 
Td to a few Kelvin). This value is larger than the torus 
radius obtained from the fi t of NIR and MIR data only 
(Ro=0.7±0.3 pc; seel lRAllI) . Finally, we can calculate the 
torus covering factor, Ct- Broader tori with more clumps 
will have larger covering factors and viceversa. According to 
our modelling, Ct = O.SioJ, which i s among the typical 
values found for Sy2 gal axies in| RAll and it is compatible 
with X-ray studies (e.g. iRicci et al.ll201in . 



4 CIRCUMNUCLEAR AND EXTENDED 
DUST PROPERTIES 

As explained above, the unresolved component of the SPIRE 
data includes the inner ring rl apart from the active nucleus 
(note the clear bump of emission at A >160 iJ.m in Figure 
[2]). To characterise the heating source of this component we 
have extracted fluxes in an aperture equals to the maximum 
size of the PSF in the SPIRE bands (i.e. the FWHM at 
500 [im; 36.9") and then subtracted the galaxy background 
emission measured in an adjacent annulus (first row in Ta- 
ble [Sjl. The blue dotted line in Figure |4] corresponds to the 
best fitting torus model (MAP), and by subtracting it from 
the latter fluxes, we can isolate the dust emission that is not 
related to the torus (red dots). The best fit that we get for 
the latter component is a greybody of emissivity e—2 and 
temperature T=28±l K, which is typical of dust heated by 
youn g stars in t i ie gal axy disk. Indeed, based on HST im- 
ages, |Buta^_et_aL| (|200J) detected ~350 diffuse bright clusters 
of Re ~11 pc, with stellar populations younger than 10 Myr. 
In a similar way, we can determine the temperature of the 
dust in the galaxy disk, by subtracting the fluxes obtained 
in the 36.9" aperture from the total galaxy fluxes reported 
in Table |3] (green squares). Figure |4] shows the fit of this 
component with a greybody of e=2 and T=19±3 K, that 
is compatible with dust heated by the interstellar radiation 
field. 



Aperture 


F70 


Floo 


F16O 


F250 


F350 


F500 


36.9" 
Total flux 


2124 
2432 


2424 
3210 


1627 
3164 


552 
1950 


187 
757 


61 
284 



Table 3. PACS & SPIRE flux densities of NGC 3081 in mjy. 
The total fluxes were obtained in apertures big enough to include 
the whole galaxy emission in each band. 
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Figure 4. Fit of the circumnuclear emission of NGC 3081 with a 
greybody of T=28±l K (red dashed line). The green dot-dashed 
line corresponds to the fit of the disk emission with a greybody 
of T=19±3 K. 



We integrated the emission of the two previous com- 
ponents and obtained their IR luminosities (1-lQQQ (im): 



Lf^'== 1.4x10'' 
equation M^ 



L© 
7.9 X 



and Lf^''=2.2xm'' 



lO-='(r/4O)-''L7fl/L0 (Mq) 



Lq. Following the 
from 



iKlaas & Elsas serl 11 19931) . w e can estimate the dust masses 
as in Radovic h et al.l 1 19991 ). We obtain a mass of 0.9x10^ 
Mq for the circumnuclear dust at T=28 K and 1.6x10^ 
Mq for the disk. The dust mass content of NGC 3081 is '-^5 
times smaller than the value reported by iRadovich et al.l 
(I999!), for the disk of the star-forming Seyfert 2 NGC 7582. 
This mass rati o is twice that of the HI masses {Mhi oc 
d^5',^(21eTw)fa: iGiovanelli fc Havnedll988r ). but it is consis- 
tent with the relative content of virial mass in both ealaxies. 
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5 CONCLUSIONS 

The FIR nuclear luminosity of NGC 3081 (on scales <1.7 
kpc in diameter) can be reproduced by warm/cold dust 
within a clumpy torus heated by the AGN. On larger 
scales (5.4 kpc), the IR emission corresponds to a cold 
dust component at T=28±l K heated by young stars 
in the galaxy disk, likely located in the rl star-forming 
ring. On the other hand, the dust located in the outer 
parts of the galaxy is heated by the interstellar radi- 
ation field (19±3 K). These components are coincident 
with the findings III studies of nearby Seyfert galaxies 
llRadovich et al.l 19991: iPerez Garcia fc Rodriguez Espinosal 



l200ll : lBendoet al.ll201ol T 



In our previous work using clumpy torus models we fit- 
ted the NIR-to-MIR SED of NGC 3081, amo ng other Seyfert 
lalajci es, with an interpolated version of the lNenkova et al.l 
20081 ) models. In this letter we have repeated the fit af- 
ter adding the nuclear Herschel fluxes to the SED. The FIR 
data provide information about cooler dust within the torus, 
resulting in a relatively large value of the torus outer radius 
and a flat radial distribution of the clumps. 
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